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ABSTRACT 

An experiment to test whether gravitation is described 
by metric or no -metric theories was carried out by B.F. 
Farrell, E. M. Ncttison, and R.F.C. Vessot of the Smith- 
sonian Astrophysical Observatory and J. Turneaure and C. 

Will of Stanford University. The frequencies of two atomic 
hydrogen masers and of three superconducting-cavity stabil- 
ized oscillators (SCSOs) were compared as the ensemble of 
oscillators was moved in the sun's gravitational field by 
the rotation and orbital motion of the earth. Metric grav- 
itation theories predict that the gravitational redshifts 
of the two types of oscillators ar identical, and that there 
should be no relative frequency shift between the oscillators; 
non-metric theories, in contrast, predict a frequency shift 
between masers and SCSOs that is proportional to the change 
in solar gravitational potential experienced by the oscil- 
lators. The results of the experiment show that 

blLL s 2.0x1 C”^ 

Aiji/c ? 

where Af/f is the relative frequency shift of masers and 
SCSOs, expressed as a fraction of their frequencies, Af> is 
the change in gravitational potential, and c is the speed 
of light; the limit represents the upper bound of the 95a 
confidence interval. The results are thus consistent with 
metric theories of gravitation at a level of 2%. 


INTRODUCTION: 


The development in recent years of extremely high stability 
frequency sources has given the physicist new tools for testing theories 
of relativity and gravitation. In the experiment described here, the 
frequencies of high stability oscillators were compared over periods of 
several weeks in order to test basic assumptions of relativity theory and 
to measure the characteristics of the oscillators themselves. The experi- 
ment resulted from a conjecture that if gravitation were described by a 
non-metric theory, clocks based upon different physical principles would 
run at different rates when the strength of the gravitational field they 
were exposed to was changed. ^ In this experiment we measured the 
relative frequencies of two hydrogen masers and three superconducting- 
cavity stabilized oscillators as the oscillators moved in the gravitational 
field of the sun as a result of the rotation and revolution of the earth. 

The output frequency of the hydrogen maser is based upon the 

hyperfine splitting of ground state atomic hydrogen. The construction and 

(2 31 

characteristics of the masers are described in detail in the literature. ’ ' 
Briefly, hydrogen atoms in the upper hyperfine level (F = 1, m p = 0) enter 
a Teflon-lined quartz sphere contained in an evacuated microwave cavity 
resonant at the hyperfine frequency (1.42 GHz) and are stimulated to 
radiate to the lower hyperfine level (F = 0, nip = 0) by the microwave field 
in the cavity. The output frequency is shifted from the hyperfine 
frequency by interactions of the atoms with the Teflon wall surface, with 
an imposed axial magnetic field, and with each other, and by variations 
in the cavity resonance frequency that can be caused by changes in 
temperature, barometric pressure, mechanical stress, and other sources. 
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The SCSO consists of an electronic self oscillator whose frequency- 

determining element is a high-Q superconducting cavity maintained at the 

temperature of liquid helium. The cavity's dimensions, and thus its 

resonance frequency (approximately 8 GHz), can be affected by such external 

effects as temperature and mechanical stress, including stress changes due 

to variations in the strength of the local gravity field and in the cavity's 

angle of tilt relative to the vertical. 

Section 2 of this report describes the theory underlying the 

experiment. The experimental procedures and methods of data collection 

are described in Section 3, and the data analysis is described in Section 4. 
The conclusions of the experiment are discussed in Section 5. 

2 . THEORY 

This experiment, which tests whether gravitation is described by a 
metric or a non-metric theory, is based upon the following conjecture: 

"Every non-metric theory of gravity predicts a gravitational redshift 

(1 

that depends on the nature of the clock whose redshift is being measured." 

(4) 

Quantitative predictions are obtained from the Lightman-Lee formalism, 
which is ipplical.e to certain classes of gravitational theories and which 
describes the behavior of electromagnetic fields and charged particles in 
the presence of an external, static, spherically symmetric gravitational 

PM 

potential U Q = — ^ • Here G is the gravitational constant, M the mass of 
rc 

the external body, and r is the distance between the body and the observa- 
tion point. The potential U 0 i s normalized to the square of the speed of 
light, c. In the limit of a weak gravitational field (U 0 <<1) the 
functions describing the particles and the electromagnetic fields can be 
expanded in terms of dimensionless parameters r , , . . . , j t^, . . . , 

whose values depend on the particular gravitational theory being used. 
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Every metric theory has r . = T. = 0 for all i; every non-metric theory 
has at least one of the r 's or T's non-zero. 

We imagine an experiment in which the frequency of a clock is 
measured as it moves from one point to another i r. a gravitational field, 
with a gravitational potential difference AU between observation points. 
The predicted redshift for a hydrogen maser clock is ^ 

( : V) c = (1-4 r + TJ AU = (1 + u m ) A U (1) 

v f ' maser v o 1 ' ' ni ' 

For a clock whose frequency is determined by the dimensions of a 
microwave cavity, such as an SCSO, the redshift prediction is 

(~T ^ SCSO ^ 1 ” ! o 2 T 1 ) A U = (1 + rt s > AU ^ 

If a maser and an SCSO together experience the potential change 
AU, their relative frequency difference will change by : 

<-/)„ - <"T-> S = ! < T 1 - 2 r o> AU = K - “s> AU = "ms AU (3) 

All metric theories (T^ = F q = o) predict that the clocks will not 
change frequency with respect to one another, that is, a = 0. 

In this experiment we attempt to place a bound on the magnitude of n . 

The largest time-dependent component of gravitational potential 
seen by the clocks is the variation resulting from the motion of the 
laboratory in the sun’s gravitational field as the earth rotates on its 
axis and travels in its orbit. The experiment was carried out in 
Palo Alto, California in late March and early April, 1978. At that time 
of year the orbital phase of the earth was approximately 90° with respect 
to perihelion, and the earth-sun separation was changing at its most 
rapid rate. Thus the potential has a component changing with a diurnal 
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period and a component varying nearly linearly with time. The expression 
for the time-varying part, U, of the potential at the position of 
Palo Alto is: 

U ( t ) = 2.84 x lO -12 (t-tj l 1 - 4.9 x 10" 5 (t-t*) 2 ] (4) 

-3-2 x 10~ 13 Cos | 2n (t-tj ] 

Here t is measured in solar days and t* = April 4.5 (that is, noon on 

p 

April 4). The term (t-t*) in square brackets results from the fact 
that the radial motion in orbit is described by a sine function with a 
period of one year; for intervals of the order of tens of days, as in 
this experiment, the departure from a linear function of time may be 
ignored. The expression for the lowest-order terms of the time varying 
gravitational potential, then, is 

U (t) = 2.84 x 10" 12 (t-tj - 3.2 x 10" 13 


Cos ( 2 it (t-tj J (5) 
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3. DESCRIPTION OF EQUI PMt'NT AND DATA ACQUISITION 
3.1. Laboratory Setup 

The experiment was performed at Stanford University, Palo Alto, 
California. Three SCSO's built at the Stanford High Energy Physics 
Laboratory were contained in a single liquid-helium dewar in a 
temperature-controlled , isolated inner room within a laboratory area. 

The dewar' s liquid nitrogen shield was cooled by a recirculating 
refrigerator using helium gas as the transfer medium. This eliminated 
the need for refilling the dewar during the experiment. The dewar was 
equipped with a servo-controlled anti -tilt device to reduce the effect 
of floor tilt on the SCSOs’ frequencies. Two model VLG-10 masers, built 
by the Smithsonian Astrophys ical Observatory (SAC), were placed in the 
outer room of the laboraco-y. Because the tanperature of this room was 
strongly affected by outside weather conditions, a window air-conciitioner 
wi t(i an on/off control was installed. To moderate the temperature 
variations caused by the cycling of the air-conditioner, the masers were 
placed in large wooden enclosures equipped with temperature controlled 
exhaust fans that helped keep the air surrounding the masers at a 
constant temperature. 

3.2 Frequency Comparison and Da t a Acqu i si ti on 

The frequencies of the oscillators were compared using the system 
shown schematically in Figure 1. (The SCSOs are referred to as SI, S2, 
and S3, and the masers as P5 and P6. Subscripts indicating frequency 
differences use numerals only: 1, 2, 3, for SCSOs and 5, 6 for masers.) 

A 5 MHz voltage-controlled crystal oscillator (VCXO) was phase-locked 
to the output of SI, and its output was multiplied and mixed with the 


* 
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signals from the four other oscillators, producing beat signals whose 
periods were measured by counters. The difference frequencies of 
S1-S2 and SI -S3 were set by means of frequency synthesizers to approxi- 
mately 5 Hz, while the S1-P5 and SI -P6 beat frequencies were set to 
approximately 1.2 Hz. The intervals over which the beat periods were 
averaged are given in Table I. 

TABLE I - AVERAGING INTERVALS FOR DIFFE R EN CE FR EQUENCY ME ASUREMENTS 


Oscil lators 

Number of Beat 
Periods Averaged 

Approximate 
Beat Frequency (Hz) 

Approximate 

Averaging Interval (Sec. ) 

S1-S2 

1500 

5 

300 

SI -S3 

1500 

5 

300 

S1-P5 

100 

1.2 

83 

S1-P6 

300 

1.2 

250 


The counters used to measure the average periods between SI and S2, 

SI and S3, and SI and P6 were zero-deadtime units that did not introduce a 
time lapse between adjacent samples. The SI -P5 counter had a one-period 
deadtime and thus introduced a 1 per-cent time lapse between 100-period 
averaging intervals. Each period measurement was recorded along with the date 
and the time at the end of the measurement. 

An additional frequency measurement system was operated simultaneously 
with the one described above. The periods of the beat signals S1-P5 and 
S1-P6 were measured by counters that measured the elapsed period for 100 
cycles, with a one-period deadtime between measurements. These data were 
recorded on magnetic tape by HP9825 calculators. Before and after the 
experimental runs, measurements over 1-, 10- and 100- period averaging intervals 
were made and recorded. In addition, before and after the runs the maser 
frequencies were compared directly by a method similar to that shown in Fig. 1, 
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in which a b MHz VCXO was phase-locked to the output of one of the 
masers. These pre-and post-run comparisons gave quick checks of the 
oscillators' stabilites and verified their proper operation. 

Measurements of ambient conditions were recorded at forty-minute 
intervals together with the date and time of each measurement. The 
types of ambient data measured are given in Table II. 

TABLE I I - AMBIENT CONDI T IONS MONJ J0RED_ DU R I N G EXPERIMENT 

SCSO CONDITIONS 

Dewar top-plate temperature 

Dewar top-plate temperature controller 

Dewar tilt, X and Y directions 

Dewar tilt-control servo currents, X and Y directions 

Dewar helium vapor pressure 

SCSO room temperature 

SCSO room temperature controller 

SCSO gunn oscillator temperatures (3) 

MASER CONDITIONS 

Maser frame temperatures (?) 

4. Data Analysis 

4.1 Data Preparation 

The data discussed here were obtained from the experiment conducted 
during April 1-11, 1978. (A preliminary run for the purpose of checking 
the experimental conditions was conducted during March 19-28. Prior to 
the April run, an anti-tilt servo-mechanism was installed on the SCSO 
dewar, and the outer laboratory air-conditioning was improved.) Data 
analyses were carried out independently at Stanford University and at SAG. 
Both analyses included examinations of the linear drift and of the diurnal 
variation of the frequency. In studying the diurnal variation, the 
Stanford effort concentrated on use of Fourier transforms of the data, 


while the SAO work used regression and correlation techniques on both 
the frequency data and the time derivative of frequency. Similar 
results were obtained with the two approaches. In this report wt 
discuss mainly the SAO data analysis. 

The four sets of beat period data were transformed Into relative 
frequency changes, fj 2 » f-|g» and fjg. (All frequency differences 
were divided by the oscillator frequencies and expressed as fractional 
differences.) Data points that were clearly Identified as being 
spurious were removed and the frequencies were Identified at Identical 
twenty-minute Intervals using a forty-minute Hanning filter. The four 
resulting time series were suitably added together to yield the ten 
frequency difference functions that are possible with five oscillators. 
These are shown In Figs. 2-11. 

The ambient condition data were Interpolated to produce time 
series whose values were in temporal registration with the frequency 
data. Of the 13 streams of ambient condition data, only three showed 
behavior that resembled variations In the frequency differences. These 
three were the vapor pressure of the helium gas In the SCSO dewar and 
the temperatures of each of the maser cabinets. They are referred to 
as HP, T5, and T6, respectively, and are shown in Figs. 12 to 14. 

During the first two days of the experiment the helium vapor pressure 
(HP) fell rapidly. The SCSO frequencies, particularly that of S2, 
showed corresponding rapid shifts. During the eleventh day of the 
experiment the liquid helium level fell below a critical level, causing 
a change in HP and an abrupt change In the SCSO frequencies. In order 
to eliminate data that were obviously affected by strong systematic 
effects, we disregarded the data from the beginning and end of the 
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experlment, keeping 8 days' of data between days 2.8 and 10.8. (Time 
s measured In ,>olar days beginning at one minute after midnight on 
me morning of April 1.) 

All sjbsequent analyses of the data use as Input these three 
ambient- condition time series and four frequency difference time series. 

4.2 Linearly Varying Component of Solar Gravitational Potential 

As shown in Eq. 5, the solar potential consists primarily of two 
terms: a diurnal variation with dally peak-to-peak change of approxi- 

mately 6 x 10" 1 J , and a linear variation with a dally change of 

-12 

approximately 3 x 10 . Because of the relatively large magnitude of 

the linear term, the existence of a gravitational effect on the 

oscillators would appear in the oscillator frequency differences primarily 

as a linear change with time. The linear variations of the frequency 

differences are examined in this section. A search for a diurnal component 

ascribable to the solar potential Is discussed in section 4.3. 

From Figs. 2 to 11 It is apparent that all of the frequency 

differences, including those for maser-maser and SCSO-SCSO pairs, have 

-13 -14 

relative drifts of up to 3.6 x 10 in 8 days, or 4.5 x 10 per day. 

-12 

Since the solar potential is varying at the rate of 2.84 x 10 per day, 
we can quickly set an upper limit on the quantity a ms of roughly 
1.6 x 10 if we ascribe all of the frequency drift to the gravitational 
potential. This estimate does not, however, take into account non- 
gravitational systematic effects acting on the oscillators. The effects 
of the helium vapor pressure and of the maser frame Temperatures were 
removed by performing a multiple function linear regression on each 
frequency-difference time series. The independent functions used in the 
regression were the calculated solar gravitational potential U and the 
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ambient effects appropriate to the particular frequency differences} 
for example, HP, Tg and U were regressed on f^; HP and U on fj 3 ; 

Tg, and Tg, and U on fgg} and so forth. Of the ten frequency functions, 
six (written f^ ) are frequency differences between dissimilar oscil- 
lators (maser-SCSO), and four (written f^) are differences between 
similar oscillators (maser-maser and SCSO-SCSO). The regression co- 
efficient, b^, o f U on each of the f^ Is proportional to any component 
of f kn that varies with the gravitational potential U, and thus gives 
a value for a ms . (Similar oscillators can have no mutual frequency 
difference Induced by U, so the regression coefficients of U on the f t 
measure the presence of unidentified systematic effects, or of ambient 
effects containing linear or diurnal components that were not completely 
removed by the regression.) The regression coefficients are glv^.. in 
Table III. 


TABLE III 


REGRESSION COEFFICIENTS OF U ON FREQUENCY DIFFERENCES BETWEEN PAIRS OF OSCILLATORS 


REGRESSION COEFFICIENT STANDARD ERROR 

OSCILLATORS |Nj ° f b ij 


5-1 

1.45 

X 

10" 2 

3.4 

X 

1 

O 

5-2 

1.25 

X 

10‘ 2 

2.4 

X 

10" 4 

5-3 

0.57 

X 

10" 2 

3.1 

X 

10" 4 

6-1 

1.50 

X 

10’ 2 

7.3 

X 

10" 4 

6-2 

1.39 

X 

10 -2 

5.9 

X 

10" 4 

6-3 

0.67 

X 

10" 2 

4.5 

X 

10" 4 

6-5 

0.30 

X 

10" 2 

1.7 

X 

10" 4 

2-1 

0.36 

X 

10' 2 

2.6 

X 

10" 4 

3-2 

0.58 

X 

10 " 2 

3.2 

X 

10" 4 

3-1 

0.93 

X 

10" 2 

3.6 

X 

10" 4 


The standard errors of the regression coefficients, which are the S' 

errors due to scatter of the data about the regression lines, are at 
most 7 per cent of the values of the coefficients, while the variations 
among the coefficients are on the order of magnitude of the coefficients; 
hence, the errors in determining the individual coefficients (the 
standard errors) do not significantly affect the variations among co- 
efficients and will be ignored in the following analysis. 

ft 

Table III shows that the regression coefficients--and thus the i 

relative drift rates--for masers verses SCSOs are all positive (that is, j 

both masers increase in frequency faster than any of the SCSOs), and that the 1 

drift rates between dissimilar oscillators are for the most part greater * 

■J 

than those between pairs of similar oscillators, indicating that the two ] 

groups of oscillators are changing frequency at significantly different rates. J 


* * * t 

... V . v ■■ ■ ■ ;• Mf. 
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The measurements of frequency differences between dissimilar 

oscillators constitute a group of experiments, each leading to a value 

for a ms • from which we wish to determlng the mean value "oT ms and Its 

confidence Interval. We cannot, however, simply use the six regression 

coefficients b.„ between dissimilar oscillators to estimate a ” 
xn ms 

because these measurements are not Independent. Since the frequency 
data are related by linear superposition, only four of the six coeffic- 
ients for dissimilar oscillators are statistically Independent. For 
example, the frequency differences fgj, fg 2 , fgj t and f fi3 produce the 
differences fg 3 and fg 2 ; by linear superposition, therefore, the 
regression coefficients bg 3 and bg 2 are related to bg^, bg 2 , bg^ and bgg. 
Of the fifteen sets of four coefficients that can be chosen from the 
six b )cn • three sets are not complete in that they do not generate the 
remaining two coefficients; these sets are (b g j , b g2 , b g1 , b g2 ) and 
(bg 2> bg 3 , b g2 , b g3 ), and (bg 3> b &1 , b g3 , b gl ). The remaining twelve 
sets of coefficients represent different ways of measuring a ms , that is, 
they can be regarded as twelve possible sets of four independent experi- 
ments. The sets are not, however, independent of one another, and we 
must choose the averaged results of one set to estimate a ms . 

Table IV gives the average regression coefficients for the twelve 

— -2 

sets of data. The maximum value of b is 1.25 x 10 and the minimum is 
_? 

0.99 x 10 , the span between these being less than the typical standard 

_2 

error of approximately 0.4 x 10 , indicating that our conclusion will 

not be sensitive to the choice of data set. 


TABLE IV 


Mean Regress -ion iflefilcients. fur. Independent Sett, ei- Jlit s.iuit lar. Qst.t 1 1 atons 


SET OSCILLATORS 


IV-1 

5-1, 

6-1, 

5-2, 

5-3 

IV-2 

5-1, 

6-1, 

5-2. 

6-3 

I V-3 

5-1, 

6-1. 

6-2. 

5-3 

IV-4 

5-1, 

6-1, 

6-2, 

6-3 

IV-5 

5-2. 

6-2, 

5-3, 

5-1 

IV-6 

5-2, 

6-2, 

5-3, 

6-1 

I V-7 

5-2, 

6-2. 

6-3, 

5-1 

IV-8 

5-2, 

6-2, 

6-3, 

6-1 

IV-9 

5-3, 

6-3, 

5-1, 

5-2 

IV-10 

5-3, 

6-3, 

5-1. 

6-2 

IV-1 1 

5-3, 

6-3, 

6-1, 

5-2 

IV-12 

5-3, 

6-3, 

6-1, 

6-2 


MEAN REGRESSION 
C OEFFI CIENT b 


1.19 

X 

10" 2 

1.22 

X 

10' 2 

1.23 

X 

10‘ 2 

1.25 

X 

10" 2 

1.1/ 

X 

IQ’ 2 

1.18 

X 

io’ 2 

1.19 

X 

CM 

1 

o 

1.20 

X 

_2 

10 L 

0.99 

X 

10‘ 2 

1.02 

X 

io -2 

1.00 

X 

10" 2 

1.03 

X 

io' 2 


STANDARD DEVIATION 


_ Of 

JL. 

.43 

X 

10‘ 2 

.38 

X 

10“ 2 

.44 

X 

10' 2 

.39 

X 

10 -2 

.41 

X 

10' 2 

.42 

X 

io' 2 

.36 

X 

io -2 

.37 

X 

io" 2 

.43 

X 

10‘ 2 

.46 

X 

10 c 

.45 

X 

IO' 2 

.48 

X 

Kf 2 


A.Jg. 


f 


**-**M*4Jt 
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To estimate an upper bound for a „ we use the maximum value for 

ms 


b and the maximum standard deviation; thus the maximum value of a 


ms 


estimated from the linear component of the solar potential is 

aJ 1nMr ■ (1-25 i' M. » - lO' 2 (6) 


ms 


The t distribution with three degrees of freedom gives the 95 per-cent 

1 inear 


confidence interval for 


°fos 


as 


0.49 x 10” 2 < ot mc linear < 2.01 x 10" 2 

— ms — 


(7) 


Hence with 95 per-cent confidence we estimate the upper bound as 


S ' ,near < 2.01 x 10’ 2 

ms — 


( 8 ) 
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4 • 3 Diurnal Component of Solar Gravitational Potential 

Because the diurnal component of U is far smaller than the 
linearly varying component, the regression calculation described in 
Section 4.2 is relatively insensitive to the diurnal component. This 
was verified by a regression calculation on the frequency data in which 
the diurnal term was treated separately from the linear term; this 
yielded no significant diurnal component of f at the 95 per-cent 
confidence level. In order to search for a gravitational effect without 
being obscured by long-term drifts, we used advanced time-series analysis 
techniques to examine the data for the presence of a diurnal variation 
in phase with the solar gravitational potential. 

The frequency differences of Figures 2 to 11 show substantial 
long-period fluctuations with durations from fractions of days to several 
days. This non -stationary behavior is characteristic of noise produced 
by a first-order auto-regressive process, that is, one with the ability 
to integrate past effects. The noise is far from white; rather, it has b' 

a distinctly red, or low frequency, character. In order to be able to 
detect shorter-period variations, we removed the low frequency noise by 
transforming the frequency functions with the inverse of the auto- 
regressive process, namely, the time derivative. Figs. 15 to 24 show 
the time derivatives of the frequency series; short-term variations with 
periods of approximately one-half day are clearly visible. (We will, 
henceforth, refer to these time derivatives of the frequency series as 
derivative series or frequency derivatives.) The auto-correlation 
functions of the frequency derivatives, as well as Fourier transforms of 
the frequency data calculated at Stanford University, showed that the 
dominant periodicity of several of the frequency series was approximately 
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0.51 days. We confirmed these period estimates with greater precision 
by the use of super-resolution, a non-linear least squares regression 
technique in which a trigometric function A cos ( ^ Is fitted 

to the data by adjusting the amplitude (A), period (T), and phase (<t>) of 
the ^unction. The period of the semi-diurnal perturbation was determined 
by this method to be 0.507 - .002 solar days. This is significantly 
longer than half a solar day and suggests the effect of the lunar 
gravitational tide. 

Possible mechanisms of lunar perturbation are changes in the local 
gravitational acceleration g (the effect that causes the ocean tides), 
and variations in the atmospheric pressure p, which is affected by the 
variation in g, resulting in the atmospheric tides. To determine the 
source of the perturbation we cross-correlated the derivatives of the 
frequency differences with the time derivative of the calculated value of 
local g and with the derivative of the barometric pressure measured at 
Palo Alto during the experiment. The closest correlation was between 
the frequency derivative and the barometric pressure derivative. The 
barometric pressure is shown in Fig. 25; its time derivative, p, together 
with that of a typical frequency derivative series, is shown in Fig. 26. 

To search for the presence of a solar diurnal component in the 
frequency derivatives, we first regressed p on the f . obtaining the 

• J 

regression coefficients of Table V. 
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TABLE V 

REGRESSION COEFFICIENTS FOR BAROMETRIC PRESSURE DERIVATIVE t> REGRESSED ON 

FREQUENCY 


OSCILLATOR 

PAIR 

1- 5 

2- 5 

3- 5 
1-6 
2-6 
3-6 
5-6 
1-2 
2-3 
1-3 


REGRESSION COEFFICIENT x 10 14 
lINCHJjgj-J 

12.0 

6.4 

- 8.1 

33.8 

25.0 

18.1 

22.6 

9.3 

7.0 

16.1 


Note that for oscillator pairs involving maser P6, the regression 
coefficients, which measure the sensitivity of the oscillators' frequen- 
cies to changes in barometric pressure, are substantially larger than 
for most other oscillator pairs. The soi.rce of this difference, and its 
possible effect on the results of the experiment, are discussed in 
Section 5.2. Using these coefficients we then removed the effect of f> 
from the f.., yielding residual frequency derivatives F. .. 0 was then 

I J * J 

regressed on the residuals to determine regression coefficients a^. 
(Because the functions in the regression were adjusted to zero mean, U 
contains only a contribution from the derivative of the diurnal component.) 
The regression coefficients a^ for dissimilar oscillators, which give 
values for ot^, are listed in Table VI. 
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TABLE VI 

REGRESSION COEFFICIENTS OF U ON RESIDUAL FREQUENCY DERIVATIVES 
( BAROMETRIC PRESSURE REMOVED ) 

OSCILLATORS REGRESSION COEFFICIENT 




a ij 

1-5 

.48 

x 10" 2 

2-5 

.59 

x 10 -2 

3-5 

.62 

x 10" 2 

1-6 

.99 

CVJ 

I 

o 

X 

2-6 

1.04 

CVJ 

1 

o 

X 

3-6 

1.37 

CVJ 

1 

o 

X 

2-1 

.04 

CM 

1 

o 

X 

3-2 

.21 

x 10‘ 2 

3-1 

.17 

x 10“ 2 

5-6 

.52 

X 

o 

1 

ro 


As discussed in Section 4.2, only certain sets of four of the six 
coefficients between dissimilar oscillators can be regarded as in- 
dependent measurements of a ms - The twelve possible sets of independent 
dissimilar oscillator pairs, and the mean regression coefficient determined 
by each set, are given in Table VII. 


X* 
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T ABLE VII 

: I ENTS a kn OLON/kn f 
DISSIMILAR OSCILLATORS 


MEAN REGRESSION COEFFICIENTS a kn OF U ON F kn FOR INDEPENDENT SETS OF 


OSCILLATORS 


COEFFICIENT a kn 


DEVIATION OF 


1 

1-5, 1-6, 

2-5, 

3-5 

.67 x 

10-2 

CM 

CM 

X 

CM ( 
1 

o 

2 

1-5, 1-6, 

2-5, 

3-6 

.86 x 

10“ 2 

.41 

X 

10- Z 

3 

1-5, 1-6, 

2-6. 

3-5 

.78 x 

10 -2 

.28 

X 

10 -2 

4 

1-5, 1-6, 

2-6, 

3-6 

.97 x 

10’ 2 

.37 

X 

10 2 

5 

2-5, 2-6, 

3-5, 

1-5 

.68 x 

10' 2 

.25 

X 

10" 2 

6 

2-5, 2-6, 

3-5, 

1-6 

.81 x 

10" 2 

.24 

X 

10" 2 

7 

2-5, 2-6, 

3-6, 

1-5 

.87 x 

10' 2 

.41 

X 

10" 2 

8 

2-5, 2-6, 

3-6, 

1-6 

1.00 x 

10" 2 

.32 

X 

10‘ 2 

9 

3-5, 3-6, 

1-5, 

2-5 

.77 x 

10 -2 

.41 

X 

10" 2 

10 

3-5, 3-6, 

1-5, 

2-6 

.88 x 

10 -2 

.41 

X 

10‘ 2 

11 

3-5, 3-6, 

1-6, 

2-5 

.89 x 

10 -2 

.37 

X 

10’ 2 

12 

3-5, 3-6, 

1-6, 

2-6 

1.00 x 

10' 2 

.31 

X 

10' 2 


The maximum value for a^ Is 1.00 x 10" z and the maximum value 
for the standard deviation of a^ is .41 x 10 ; thus the maximum 

value of "aj J)S estimated from the diurnal component of the solar 
potential is 

~“ms diurnal = (1.00 t .41) x 10 -2 (9) 


The 95 per-cent confidence interval for diurnal 


35 x 10 


< T diurnal < 1.65 x 10' 2 
— ms — 


and the estimated upper limit for a. diurnal i s 

r ms 


< 1.65 x 10 


diurnal 


(ID 
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5. DISCUSSION OF RESULTS 
5.1 Linear Drift of Frequen cy 

We have seen that the masers have frequency drift rates of several 
14 

parts in 10 per day relative to the SCSOs, with the maser frequencies 

consistently increasing with respect to the SCSOs. If all of this drift 

is ascribed to the solar gravitational potential one concludes, as in 

-2 

Section 4.1, that a has a value of 0 to 2)x 10 . However, frequency 

standards are often subject to long-term systematic drifts; in particular, 

it has recently been found at SAO and at other laboratories that 

some hydrogen masers have slow upward drifts in frequency at a rate of 

14 

up to 1 to 2 parts in 10 over time spans of many months. This is 
believed to be due to settling of the ground joints of their Cervit^ 
resonant cavities. Over the year prior to this experiment, P6's resonant 
cavity decreased in length at a rate corresponding to a change in output 
frequency of 1.6 x 10”^/day. (This change was compensated by retuning 
the maser periodically, that is, adjusting the cavity resonance frequency 

(21 

so that the output frequency is independent of the hydrogen beam intensity. 

This procedure returned the output frequency to its canonical value; 

between tuning, however, the frequency increased approximately linearly.) 

Data on maser P5's cavity frequency drift are fragmentary; however. 

Table III shows that the 5-6 regression coefficient, and hence the drift rate 

between masers P5 and P6, is small compared to those of other oscillator 

pairs, leading to the conclusion that P5's frequency had a long-term 

behavior similar to that of P6. It appears, then, that the frequencies 

-14 

of both P5 and P6 drifted upward at rates of approximately 1 to 2 x 10 / day 

during the experiment due to resonant cavity deformation. A spontaneous 
frequency drift of 1.6 x 10 ^/day (P6's drift rate) due to the masers 
contributes to « ms approximately 0.54 x 10 , 


which Is 40 percent of the value of a ms obtained In Eq. 6. 

o linear f 0r this cavity- Induced drift gives 
ms 

- linear ( corTec ted) ■ (0.71 - 0.48) x 10 * 
ms 

and the upper bound for the 95 percent confidence Interval 


Correcting 


( 12 ) 


becomes 


a 


ms 


Hnwr (corrected) < 1 .47 x 10' 2 . 


03) 
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5.2 Diurnal Frequency Variations 

As shown by the regression coefficients of Table V, changes 

In barometric pressure affect all of the frequency differences, 

particularly those Involving SCSO SI or maser P6. This Is caused 

by the physical construction of the oscillators. The SCSOs' 

barometric sen:1t1v1t1es arise from the fact that they are suspended 

by waveguides from a plate covering their dewar; changes In barometric 

pressure deflect the plate, which tilts the SCSO cavities, causing 

changes In cavity dimensions that alter the cavity resonance 
(91 

frequencies. ' ' 

The masers' barometric sensitivities were measured In a pressure 

chamber at the Jet Propulsion Laboratory aft ir this experiment. P6's 

- 1 3 

sensitivity was found to be approximately 3x10 /Inch Hg, which is 
consistent with the results shown In Table V and Is far greater than 
the sensitivities of other masers. When P6 was subsequently d 1 s- y 

assembled, the anamalous barometric sensitivity was found to have been 
due to a mechanical interference within the maser, which was then 
corrected. In contrast, P5 has had no detectable barometric sensi- 
tivity In tests at SAO; measurement uncertainty places an upper bound 

- 1 5 

of approximately 8x10 /Inch Hg on P5's sensitivity. 

The oscillators' barometric sensitivities can affect the estima- 

$ 

tion of the diurnal solar effect in the following way. The local varia- , 

tion of g, and thus the atmospheric tides, have a diurnal component 
in addition to the dominant semidiurnal term. Because the 

regression analysis does not yield a perfect fit of the barometric 
pressure derivative p to the frequency derivatives f^, some of the 
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dlurnal variation due to the pressure remains In the frequency residuals 
Fj j after removal of p. Residuals with large remaining diurnal pressure 
variations may then have high correlations with U, leading to high 
estimations for o ms . This Is suggested by Table VI: the regression 

coefficients a k g (k ■ 1,2,3) for P6-SCS0 pairs are approximately twice 
as large as the a kg . However, the regression coefficients for similar 
oscillators contradict the notion that residual pressure-induced diurnal 
variations are the major source of correlation between the frequency 
residuals and the solar gravitational potential. Although the barometric 
sensitivities for SCSO-SCSO pairs are at least as large as those for 
P5-SCS0 frequency differences, the a rt for SCSO-SCSO frequency residuals 
are less than half the magnitudes of the a^g. Also, agg has the same 
magnitude as the a^g even though fgg has a much larger sensitivity than 
do P5-SCS0 frequency differences. 

This behavior can be explained by the effect of the maser room 
air temperature, which, as shown In Figs. 13 and 14, has strong diurnal 
variations. These variations are due to the dally changes of outside 
air temperature, which are imperfectly compensated for by the laboratory 
temperature control system. (The SCSOs were held at a constant tempera- 
ture by the dewar and by careful temperature control cf the Inner 
laboratory room.) The effect of maser temperature was not removed In the 
analysis of the diurnal variation, primarily because the temperature 
functions are not orthogonal to the solar potential and thus could not 
be removed unambiguously by a regression technique. The mechanical 
Interference in maser P6 that caused Its anomalously large barometric 
sensitivity was capable of producing a large frequency sensitivity to 
ambient temperature, resulting In large correlations between U and the 


-24- 


frequency residuals Involving P6. Because P5 and P6 experienced 
almost Identical ambient temperature variations, their frequency 
difference would be less affected by temperature than would P6-SCS0 
frequency differences. This would account for the smaller value of 
agg compared to a^g. The effect of temperature could properly be 
removed by measuring the masers' temperature sensitivities In a 
thermally-controlled chamber to determine the transfer function 
between ambient temperature and maser frequency. By applying this 
transfer function to the maser temperatures measured during the ex- 
periment, the temperature effect could be removed as a calibration, 
rather than by means of an adjustable parameter regression. This 
technique would minimize the 1 1 kl 1 hood of confusing the effect of 
temperature with a diurnal frequency variation due to the solar 
gravitational potential. 

CONCLUSION s 

The results of this experiment are consistent with the predictions 
of metric theories of gravitation, with the maximum deviation of a m$ 
from zero being 

a mc < 2.01 x 10" 2 
ms — 

This limit is the largest value of the upper bound of the 95 percent 

confidence interval estimated from the frequency data uncorrected for ^ 

long-term frequency drifts. After correction for cavity-induced maser 

frequency drift, the frequency data give a m$ ^ near (corrected) 

_2 

<_ 1.47 x 10 , while the diurnal component with barometric pressure 
sensitivity removed gives a ms ^^ urna ^ <_ 1.65 x 10~ 2 . 
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FIGURE CAPTIONS 


1. Block diagram of oscillator frequency comparison system. 

2-11. Relative frequency differences between pairs of oscillators 
during experiment. 

12. Helium vapor pressure In SCSO dewar during experiment. 

13. Temperature of outer frame of maser P5 during experiment. 

14. Temperature of outer frame of maser P6 during experiment. 

15-24. Time derivatives of frequency differences between pairs of 

oscillators during experiment. 

25. Barometric pressure at Palo Alto, Ca. during experiment. 

26. Time derivative of barometric pressure and time derivative of 
frequency difference between maser P6 and SCSO SI during experiment. 
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